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ABSTRACT 



We present the first calculation of the ^^-Wigner-Kirkwood corrections to a rel- 
ativistic system of fermions in the presence of external scalar and vector potentials. 
The method we propose allows to compute efficiently semiclassical corrections to 
one body operators such as mean energies and local densities. It also preserves 
gauge invariance and produces explicitly convergent results despite some apparent 
divergencies at the classical turning points. As a byproduct we obtain the h'^ cor- 
rections stemming from the polarization of the Dirac sea. We apply our results to 
the relativistic a-u Lagrangian in the Hartree valence approximation. We compare 
the semiclassical expansion with the exact result and with the Strutinsky average 
whenever it can be obtained. We find that the corrections are much smaller than 
typical shell effects. Our results provide convincing arguments to neglect higher 
than second order h effects in the Wigner-Kirkwood scheme to relativistic nuclear 
physics in the mean field approximation. 
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1 Introduction 



Relativistic mean field theories for nuclear systems have already a long history 
and have found a wealth of applications [0, ^. They use as a starting point a 
relativistic quantum field Lagrangian which hopefully includes the relevant degrees 
of freedom, i.e. nucleons and several mesons. In the Hartree approximation it is 
possible to reproduce the saturation properties for nuclear matter and the shell 
structure of single particle states is accurately described due to the correct sign 
and size of the spin-orbit interaction. Many of the calculations are performed in 
practice for spherical closed shell nuclei, where shell effects are expected to be least 
significant. For these systems many authors have proposed to employ a relativistic 
Thomas-Fermi (TF) model § and generalizations thereof §, § g |, g, |TU|, |n|, 0, |T3|. 

From a purely computational point of view, the study of the semiclassical expan- 
sion for relativistic systems is far from being a trivial generalization of well-known 
non-relativistic methods which have been discussed at length for nuclear systems in 
the literature (for reviews see e.g. |14, 15| and references therein). The main problem 
is that whereas in the Schrodinger equation the Planck's constant h appears in the 
second power, in the Dirac equation h appears in the first power. The final results 
for bulk properties do not depend, however, on the sign of h due to parity invariance. 
This simple circumstance makes the direct generalization of non-relativistic methods 
such as e.g. Wigner transformation or others quite time consuming and generate 
extremely complicated expressions in intermediate stages of the calculation. An ad- 
ditional difficulty in the relativistic case is the appearance of four dimensional spinor 
matrices. One should also say that the treatment of nuclei requires at least a scalar 
potential in addition to the usual electrostatic potential of atomic physics, therefore 
making the computations much more involved. This probably explains why until 
now only the corrections up to fi^ have been computed in relativistic nuclear systems 
P, B] and up to h'^ in atomic systems |TB[ . In this respect one should mention that 



in the non-relativistic case the and have been known since a long time |T 



and quite recently h -order corrections have been determined |T^. The question 
arises whether there are reasons to neglect higher orders than fi^, besides the desire 
to keep the problem at a manageable level. 

The answer to this question cannot be given in a closed form, because in general it 
is not known whether the semiclassical series is convergent or only asymptotic. In the 
absence of exact theorems only numerical investigation remains. As a practical rule 
there are cases like e.g. the non-relativistic harmonic oscillator where the Strutinsky 
average is known to reproduce the complete Wigner-Kirkwood (WK) result . 
This rule seems also to be fulfilled for the relativistic harmonic potential [|lT| in order 
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WKft^.Q However, in the cases of relevance for nuclear physics, i.e. short range and 
finite depth potentials, the Strutinsky average is in general not well defined, although 
nothing prevents from an order by order semiclassical calculation. 

In the present paper we propose an efficient computational method which allows 
to compute up to WK/i^-corrections of the semiclassical expansion for one-body 
operators following the method proposed in previous work |T^. This method is 
fully equivalent to the standard WK expansion |21|, |22|, ^ for one body operators; 
the main advantage is that one does not need to evaluate non-diagonal elements 
of the Green's functions and traces back all physical quantities to the cumulative 
number of states. As an aside we compute analytically the WK^^-corrections coming 
from the polarization of the Dirac sea. We apply our results to the relativistic a-uj 
model in the Hartree valence approximation and make a numerical estimate of the 
WK^ -corrections to the fermionic contribution to the total energy. To set in the 
proper scales we adjust the parameters to reproduce the nuclear matter saturation 
properties and the charge mean squared radius of ^°Ca P^. In doing so we do 



not pretend to make a full account of finite nuclei properties, since there are more 
realistic generalizations of this model to include Coulomb repulsion, tensor force, 
etc. [0, ^ |[. Nevertheless we feel that the general features obtained in the present 
study should not be substantially modified by these extensions of the a-uj model. 

The outline of the paper is as follows. In section 2 we review the properties of 
a D-dimensional non-relativistic case to illustrate the basic mathematical objects 
such as the determinant of the Schrodinger operator and our method to perform 
semiclassical expansions in the simplest situation. As an example we apply our 
formalism to the three dimensional harmonic oscillator. We also discuss how a re- 
ordering of the semiclassical series proposed recently p5|, ^ can be incorporated 



into the formalism, proving extremely convenient for numerical computations. Sec- 
tion 3 studies the relativistic case and the particular features of the determinant of 
the Dirac equation and explicit calculations up to ^^-order are also presented. As 
an illustration we apply and discuss those corrections in the analytically solvable ex- 
ample of relativistic harmonic potentials with special emphasis on the distributional 
character of the quantities involved. We also show an explicit calculation of the 
(finite) WKfi^-corrections arising from the polarization of the Dirac sea. In section 
4 we quickly review the a-uj model together with the parameter fixing and give nu- 
merical estimates to the WK^^-corrections to the Hartree energy. We also provide 
a prescription how to avoid the distributions in practical calculations. Finally, in 
section 5 we draw our conclusions and present some perspectives for future research. 



^The notation emphasizes that we are performing a WK expansion throughout, rather than an 
extended TF, that is, we expand in gradients of the potentials instead of gradients of the densities. 
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2 The Number of States Formalism 



2.1 Cumulative Number of States 

In this work we will consider objects which involve single sums, i.e. one body 

operators. All the information concerning them can be deduced from the single 
particle level density. As it will become clear below, it is more convenient to deal 
with the cumulative number of states. Let us consider the eigenvalue problem 

H^na{x) = En1pn^{x) (1) 

where H is the single particle Hamiltonian, are the energy eigenvalues, pos- 
sibly Qn-iold degenerated, and il^nai^) are the corresponding cigcnfunctions (a — 
1, . . . ,gn)- To fix ideas we will assume the Hamiltonian to be of the Schrodinger 
form 

although many of the features presented here can be generalized. Here V represents 
a covariant space derivative = Vj + jrAi with Ai the vector gauge potential which 
includes the coupling constant. The cumulative number of states N{E), counting 
the number of states below a certain energy E, is defined as 

oo 

N{E) ^J29nQiE-E^)^Tre{E-H), (3) 

n=0 

where stands for the possible degeneracy of the system. From here the level 
density can be derived to be 

9{E) = —N{E) = j:g^6{E-E^). (4) 

«^ n=0 

Both N(E) and g{E) are to be interpreted as distributions. For our purposes it is 
extremely convenient to consider the following representation of the step function 
distribution 



Q{x) — lim — Imlog(— X + i e) 

£— »0+ TT 



1 if x > 

I ifa; = (5) 

if X < 



the limit e — 0"*" being always understood at the end of the calculations. The branch 
cut for the log function is considered to be along the positive real axis as follows 

log2; = log + i arg2;, 0<arg2;<27r, (6) 
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which has the property 



^ log [{-X + ie)'] = log{-x + ie) . (7) 

This property will be used later in the relativistic case. With these definitions the 
cumulative number of states can be written as 

N{E)^-lmTr\og{H-E + ie). (8) 

TT 

If we further use the identity 

TrlogA = logdetA (9) 
with A a linear operator, the following formula is obtained 

N{E) ^-argdetiH - E + ie), (10) 

TT 

i.e. the cumulative number of states is proportional to the phase of the determinant 
of the operator H — E + i e. It should be mentioned that the definition of N{E) 
implicitly assumes that the Hamiltonian is bounded from below. If this is not 
the case one can still proceed with the aforementioned formulas but subtracting a 
suitable infinite but energy independent constant. This point will be discussed when 
dealing with the Dirac equation. 



2.2 Physical Quantities 

Any single sum can be obtained directly from N{E). For instance, the total number 
of particles below a chemical potential fi is given by N{E — //). Similarly, the total 
energy of a system of fermions can be expressed as 

oo „ 

^r(/i) = E^n^ne(/i-K)= dEg{E)EQ{fi-E). (11) 

n=0 •' 

In this particular case integration by parts yields the simpler result 

Et{i^) ^ liN(ii) - J dEN{E)e{n- E), (12) 

which is particularly interesting when evaluating the free energy. One body local 
quantities, such as particle densities and energy densities can be obtained from suit- 
able functional derivatives with respect to external potentials with proper quantum 
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numbers. In the non-relativistic caseQ the particle density of a system of N[n) 
fermions is 



P{X^ f^) = J2J2 ^*na{x) ^na{x) &{n - En). (13) 

n=0 a=l 

Considering the variation of N{E) in eq. (^) under a variation 6V{x) of the potential 

6N{E) = - Tr {6{E - H)6V{x)) (14) 
one readily obtains the following useful result^ 

= - E E rnai^) ^na{x) 6{E - K), (15) 

yielding 

p(x,/i) = -/rfE^^e(/i-E). (16) 

The total energy density 

OO gn 

Pe{x,^i) = E E i^*na{x) Eni)na{x) Q{^1 - En), (17) 
n=0 a=l 



can be obtained in a similar way, 



PE{x,fi) = -J dE^-^^EQifi-E). (18) 

For the kinetic energy density (total minus potential) one has 

PKix, ^) = -jdE {E - V{x)) e(/x - E). (19) 

Other local densities can be also computed by introducing an external potential with 
suitable quantum numbers, for instance the current density could be evaluated by 
a functional derivative with respect to the external vector potential. 



^ The results for relativistic systems will be presented in section |[ 

^Note that the result is base independent due to the sum in the eigensubspace. 
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2.3 Semiclassical Expansion in the Non-Relativistic case. 

As we have said any one body operator can be deduced from the cumulative number 
of states which in our regularization turns out to be proportional to the phase of the 
determinant of H — E + ie. Our aim is to perform a semiclassical expansion of those 
objects. Besides an overall factor this expansion correspond to a derivative expansion 
in the potential V{x). Therefore in most of the expressions the ^-dependence can 
be directly read off. To perform this expansion we use the method of Chan ^ 
recently extended by two of us and also considered in previous work in 
the present context. In the appendix we review for completeness the main features 
of the method and provide some results which will be extensively used along this 
paper. 

Many of the features which occur in higher order WK^ corrections within our 
formalism are best exemplified in the non-relativistic case. Thus, let us consider 
the semiclassical expansion of the cumulative number of states for a system of non- 
relativistic fermions with spinor degeneracy u interacting with a spin independent 
mean field potential (the generalization to spin dependent forces is straightforward 
although somewhat more involved). This case is well known from the study of non- 
relativistic atomic systems in the context of density functional theory |^ and we 
only show it here for pedagogical reasons. We will do this up to fourth order in h 
in D-dimensions. This particular example is also interesting because it illustrates, 
in the case D = 3, the meaning of apparently divergent integrals. This feature does 
not show up at orders in h strictly lower than four. As shown in the appendix it is 
convenient to consider the function A given by 

A(x,k) = — — ^ . (20) 

Furthermore we use the notation (see caption of table ^ 

Ai = hViA 

Fij = hiViAj-VjAi) 

F^ = Fi,F,, 1,3 = 1,2,..., D, (21) 

where Fij represents the field strength tensor of the vector gauge potential Ai times 
%. In this particular case A possesses a trivial matrix structure, which simplifies a 
lot the actual calculation. One has after application of table ^ the result 

N{E) = No{E) + N2{E)+N^{E) + ... , (22) 

where 

No{E) = i Im / log(A;2 + 2m(\/(x) - E) + ze) 
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1 

N2{E) = -Trim 



N,{E) = — 



2r 



TT D(D + 2) 



Im 



(VA)2 



-(VA)^+ 



+ 2(A(V'A)f-:^A^ 



(23) 



It should be noted that the N2n{E) term contains the power ^ . To proceed 
further we make use of the integrals 



- Im / d^kk' log(A;2 + 2m{y{x) - E) + le) 

TT J 



„ D 

2ttt 



Im 



TT 



d^k 
(2^ 



1 r(. + f 



^ (2m(E-y))^ e(E-y) 



(2m) ^ 



X 



n 



n— s— 1 



{E-vp-^e{E-v) 



(24) 



yielding the following expression for the different contributions to the cumulative 
number of states 



A^o(^) 

A^2(^) 

N^{E) 



i 



2m\ ? 



J d^x{E-V)^e{E-V) 



(4^)fr(f + 1) 
(47r)fr(f) yn^J 

(47r)fr(f) U 
1 
2 



rf"^^^5|[(^-^)^-'e(E-v^) 



+ i(vV)2a|-^F2a^J (E-v)T-ie(E-y) 



One can reduce the number of energy derivatives by using the identity 

{V,V) ±f(E-V)^-VJ{E-V). 



(25) 



(26) 
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Integration by parts yields the result 



iVo(^) 
N2{E) 
N,{E) 



V 



/2m\ 2 
(4vr)#r(f + l) \ W 



V f2m\ 2 



(V2\/) 



12 



(E-i/)T-ie(E-i/) 



(1m\ 2 



(4vr)fr(f) V;. 



o 94 

+ - (V.V,\/)2 d\ - -^F' Oe} {e- v)^-' e{E - V) 



(27) 



It is interesting to notice that these formulas produce finite results if the deriva- 
tives with respect to the energy are evaluated after the space integrals have been 
performed. If the opposite process were considered, explicit divergencies would ap- 
pear due to the singularity caused by the classical turning hypersurface. After careful 
regularization the result would coincide with that obtained by using the previous 
prescription. This point can be clearly understood when dealing with a particular 
example (see next section). 

According to the general formula ([121) one can use the expressions (^) to perform 
a semiclassical expansion of the total energy of N{fi) fermions 



Et(/x) = Eo(/i)+^2(/i)+^4(/i) + ... 



(2^ 



where 
Eo(/x) 

E2{fi) 



/iA^2(/i) 
V 



(|?)"/rf^x(^-V)^+^0(;.-r) 



(4vr)Tr(f + 2) 



/2m\ 



E^{[i) 



(4.)#r(f) \ 



(4.)* r (f ) 



12 

2m\ 2 ^ 



(/i-V)T-l0(;.-r) 
1 



J J 288 l-^ ' ' 



24 



(vv) d, + - iv.v.vyd' - -,F' (^ - vp-' e(M - v) 



(29) 
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Similarly to the cumulative number of states one has that E2n{lj) contains a power 



2.4 Application to the Non-Relativistic Harmonic Potential 

To illustrate certain aspects of the fourth order semiclassical expansion not present 
in lower orders we consider the three-dimensional non-relativistic isotropic oscillator 
defined by the central potential 

V{r) = ^muj^r^ (30) 

and Ai = 0, with eigenenergies and degeneracy given by 

/ 3\ 1 

En = hiu i^n + -j , gn = -in + l){n + 2), (31) 

respectively. Hence the cumulative number of states is given by 

N{E) = uJ2-{n+l){n + 2)Q{E - E^) (32) 

n=0 ^ 

where z/ represents the spin degeneracy factor. The actual calculation of successive 
h corrections to N{E) is simplified by direct application of formula (0). This yields 

N,{E) = J^Awi(B). (33) 

Notice that the result is finite but N{E) has to be understood as a distribution in 
E. To understand how these expressions arise let us consider the following integral 
appearing in the fourth order contribution 



1/7/ ^OO 

Ni(E) = — — hudl / dxxWE -x^e(E -X 
2407r JO 



2407r 
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If we had derived with respect to the energy inside the integral, then exphcit diver- 
gencies would have occurred. After careful regularization the sum of the divergent 
terms would have produced the previous well defined result. The necessity of tak- 
ing the derivative outside the integral is not directly linked to the appearance of 
a genuine distribution in the final formula. For instance, if only two derivatives 
with respect to the energy were involved in the previous calculation, the final result 
(proportional to Q{E)) could only be obtained by our prescription or suitable reg- 
ularization. Notice that if the step function is ignored when deriving the integrand, 
then the final result is explicitly divergent. 

Similarly one can also compute the semiclassical expansion for the total energy 
of a system of N{n) fermions as defined by (|1^) obtaining 



8 \hujj 
1 7 

Mf^) = (35) 
2.5 Reordering of the Semiclassical Expansion 

From the knowledge of the number of fermions N{fi) below a chemical potential n 
and the total energy E{fi) one can obtain the energy of a given number of particles 
A, by properly adjusting /i so that A = N{fi). In what follows we will discuss a 
procedure within the semiclassical approach, which is advantageous for practical 
calculations [^, This method can be applied both in the non-relativistic and 



the relativistic cases. Our starting point is the former semiclassical expansion for 
the number of particles and the total energy 

N{fi) = No{^I) + N2{fi) + N^ifi) + . . . 

E{fi) = Eoifi) + E2ifi) + E,{fi) + . . . . (36) 
From the normalization condition A = N{fi) it is clear that fi inherits a certain h 



dependence. We proceed along the hues proposed in previous works |]26|, |2^, by 
further expanding the chemical potential into powers of h 

jj = + H2 + ^4 + ■ ■ ■ , (37) 
with the lowest order condition 

No{iJ,o)=A. (38) 
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It should be clear that by doing so we are simply reordering the semiclassical series. If 
the expansion is truncated at a certain finite order in the result of this reordering 
is consistent up to higher h corrections. One should also mention that for this 
expansion to be a true ^-expansion, i.e. that /io is independent of one has to 
declare that A is of order since this is the power which appears in Nq. We will 
always assume this particular order counting. If we insert equations and ( |5BD 
back into the first line of equation (|36|) the successive terms in the expansion of /x 
are automatically fixed giving 

Ho : Nq{iiq) = A 



Similarly one can expand the total energy to give 

Erifi) = So{fio)+S2{fio)+SM + O{h'^~''), (40) 

where 

£2ifJ'0) = ^2(/io) +/i2(/io)^o(yUo) 



(41) 



If we now apply the relationship between the number of particles and the total 
energy as given by eq. ( |l2|) one obtains after deriving with respect to fi 

Using this result we get the final simple expressions 

1 m^^o)? 



(43) 
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The advantages of this reordering in the total energy are manifold. First, the chem- 
ical potential has to be adjusted only once, namely at lowest and simplest order. 
Second, there is a clear order by order separation, that allows to treat the corrections 
in the energy additively. Finally, big cancellations are avoided in practical calcula- 
tions since the cumulative number turns out to be extremely sensitive to variations 
in the chemical potential. 

3 Relativistic Case 

3.1 Cumulative Number for the Dirac equation 

In this section we will consider the semiclassical expansion of the Dirac Hamiltonian 
with external vector and scalar fields. The results derived here will be applied to 
estimate higher WKh corrections in a mean field relativistic nuclear model. 

The Dirac eigenvalue problem of a particle in D dimensions can be written as 

H = a-P + P^x) + V{x), (44) 

where a ■ P = J2iLi a^P* with a* = 7^7* and /9 = 7°. The gamma matrices satisfy 
the usual anticommutation relation {'jfj,,ju} = '^g^v where the Minkowskian metric 
is (1,— 1,...,— 1) and the generalized momentum is given by 

P' = -ihdi- A\x). (45) 

The time independent potentials $(a;) and V{x) represent the scalar field and the 
time component of a (D-l-l)-vector field respectively. The rest mass of the particle 
m is included in the scalar field through the asymptotic condition $(00) = m. The 
cumulative number of states is defined as 

7V^(E) = Y.gn Q{E - En) = Tr e(^ - H), (46) 

n 

where Qn represents the possible degeneracy of the eigenvalue En- The Dirac Hamil- 
tonian is not bounded from below and hence the cumulative number becomes in- 
finite. Nevertheless one can redefine this quantity by suitably choosing the energy 
for which the number of states is counted from. This operation corresponds to a 
certain normal ordering. In the free case, the Dirac equation exhibits an energy gap 
of width 2m, which distinguishes between positive and negative energy eigenvalues. 
In the interacting case it seems reasonable to assume that, if the fields V and $ are 
not too strong compared to the rest mass m, some gap in the Dirac spectrum exists. 
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For the potentials which usually appear in relativistic nuclear physics this turns out 
to be the case. We will see that in the semiclassical limit there might appear an 
energy gap if the interval [max{K(x) — min{\^(x) + is not empty. The 

former discussion suggests the subtraction point to lie within the gap 

N{E) = iV^(E) - iV^(^gap). (47) 

Furthermore, it is convenient to split N{E) into the regions above {^+) and below 
the gap as follows 

N{E) = N+{E)-N4E) 

N^{E) = Yl 9rME-Er,) 
n e r2+ 

N_{E) = Yl 9nQ{En-E). (48) 

n e r2_ 

It should be noticed that with this definition both N^{E) and N_[E) are non- 
negative numbers. The level density is then given by 

'^^^ = -dE- = -lE- 

= Y.9n S{E - E^) , (49) 

neft 

where the irrelevant constant N^{Egs,p) disappears. To evaluate the semiclassical 
expansion of N{E) we will deal directly with N^{E) and then we will proceed to 
perform the subtraction and the separation between the two branches (above and 
below the gap). 

3.2 Independent Particle Approximation 

Notice that all manipulations above are intrinsic to the Dirac spectrum and do 
not yet define a many body picture. In the spirit of the Hartree approximation 
(independent particle model) we introduce a chemical potential as we did in the non- 
relativistic case. However, the total energy computed as the sum of all eigenstates 
below the chemical potential is an ultraviolet divergent quantity, which cannot be 
made finite by simply subtracting the vacuum contribution from the Dirac sea. In a 
renormalizable field theory, one could subtract additional counterterms to the total 
energy. To keep the line of reasoning straight we relegate the discussion of the Dirac 



sea corrections to section 3.5. For most of this paper we will work in the Hartree 



valence approximation, where many calculations have been done. As we will show 
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the semiclassical Dirac sea corrections can be formally obtained from the valence 
contribution. From the point of view of the many body theory this approximation 
corresponds to deal with a non-relativistic many body problem using relativistic 
kinematics. 

The valence energy of a system of A-fermions with energies between the chemical 
potential and the energy gap is given by 

EriiJ,) = 1^ N{fi) - JdE N{E) e(/i - (50) 

with the normalization condition 

A = N{i2). (51) 

In evaluating several single particle densities, it is necessary to make use of a gen- 
eralization of the eq. (|Tip to the relativistic case 

= T E t^Lix)^nc.{x)6{E-E^) (52) 

O^yX) ngQ±a=l 

-jirV- = T E T.^Ux)f3^Pna{x)6{E-E^). (53) 

Thus, the valence energy density is given by 
pE{x,fi) = p^(x,/i) -p^(x,/i) 

p|(x, /i) = EE ^L(a;) En 1pna{x) (± (/i - En)) 
n e U± a=l 

= TjdE^-^^Ee{±{f,-E)), (54) 

the valence particle density is 

p{x,fi) = p+(x,/i) - p"(x,/i) 

p±(x,p) = E T.^iM^na{x)e{±{fI-En)) 
n e ft± a=l 

= ^h^ww""^^^'-^^^' '''' 

and the valence scalar density reads 

ps{x,fi) = pt{x,fi) - pj{x,fi) 
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n G Q± a=l 



(56) 



The kinetic energy density could be readily obtained from the former ones by means 
of the formula 

= + (m - $(x))p+(x,;u) - (m + r(x))p+(x,;u). (57) 



3.3 Semiclassical Expansion 



Following previous work |T2| we can apply a semiclassical expansion to the relativis- 
tic case. A direct application of the logarithmic regularization to the step function 
gives 

1 



Im Tr log(if — E + ie) 



— arg det{H — E + ie). 



IT 



(5^ 



which does not allow yet to employ the general method described in the appendix. 
To do so, we consider the transformation of the Dirac operator in an odd number 
of dimensions 

= -a- P - P^x) + V{x) , (59) 

where 75 = — 'i7o7i---7D = 75^^- Using the properties of the determinant one 
obtains that 



whence 



argdet(if - E + ze) = - argdet [{H - E + ie) {H - E + ie) 



N^{E) = — ImTrlog(-p2 - U{X)) 

2lT 



(60) 
(61) 



where the following definitions have been introduced 

U{X) = <l>^-{E-V-ie)^ + iha-V{V-P<^)-^aijF,j 



F 



(62) 
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The former equation is now ready for a semiclassical treatment, since it involves 
an elliptic operator with some spinorial structure. Up to fourth order in the WK 
expansion a direct application of table |^ yields 



D{D + 2) 
1 

h 



2 A V'A + ((V.A) (V,A))^ - 2 (VA)^ 



- ^ {F,,A^y - ^ iF,^ A(V,A) (V, A) A 



(63) 



where the matrix propagator A(a;, k) ^ = k'^ + U{x) has been defined. This operator 
contains some h dependence, so that we reexpand around 

Ao{x,k)-^ = k^ + <^^ - {E -V -tef. (64) 

Direct calculation of the Dirac traces together with integration by parts yields 



iv^(E) = N^{E) + n!^{e) + N^{E) + e(;i^-^) 



(65) 



where 



<(E) 
N^{E) 



Im 



T{D) 

27r 
T{D)h^ 



27r 



Im 



d^xd'^k 

d^xd^k (1 r 



k'^ 

+ ^ (VA. 



2 I 
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(66) 
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Here T{D) is the spinor space dimension {T{1) — 2; T{3) = 4). To work out the 
momentum integrals we note that any power of Aq can be obtained as 



A" 



(-1) 



n— 1 



r{n) 



(67) 



We also use the identity 



lim — Im \og{x — iye) — 1 + sign(|/) Q{x) 

e— ►0+ TT 

and define the local momentum 

p'^{x) = {E-V{x)f - ^xf. 



(68) 



(69) 



It is interesting to notice that within a semiclassical treatment one can always define 
a gap if the interval [ma.x{V (x) — ^{x)} , mm{V{x) + ^{x)}] is not empty. In practice 
this allows to explicitly distinguish the contribution of states below and above the 
gap in the semiclassical limit by means of the formula 

1 + sign(£; - V) e(/ - k"^) = 
= 1 + e(E -V - + A;2) -Q(-[E-V)- + . 

The momentum integrals which appear in the semiclassical expansion have the form 
1 



Im 



TT 



,n-l 



(47r)^/2r(n)r 



an 
0^2 



X 



1 + Q{E -V - + k^) -q{-[E-V)- V^^TP) 



If we further subtract the contribution from the gap and explicitly separate the two 
branches we get for n > 



where 



(-i)"-^rfs + §) 



(47r)^/2r(n)r(f 



p 



e{±{E - y) - $) 



(70) 



(71) 
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The zeroth order term can be evaluated following similar steps. After suitable 
integration by parts, the final result for the cumulative number of states reads 



nHe) 



T{D) 



2 {AT^h^)DI^V{^ + 1) 



T{D) 



(v$) + {E-v) (vy)]2a|.}p^-2e(±(E - F) - $) 



3 

T{D) 



2 {A7Tn^)D/^T{D/2) 



-[$ (V*) + {E-V) {VV)]'dl + -[(V2$)2 - (VV)^]9|. + 



{vvy + (v$)^ - 6 {vvyiv^y + 4{vv- v^y ai^ + 



1 

+ 24 

+ - (V$)'][(V$)2 + $ (V2$) - (VF)2 + - y) (VV)]a|2 + 



+ 



6 



(V'<l>) (V,$) - (VV) {V,V) [<I>(V,$) + - V){V,V)] dl^ + 



+ l[$(V$) + (E-\/) (VF)f X 

X [(V$)' + $ (V'$) - (WVy + iE-V) (VV)] ^2 + 

+ ^ [(V$)' + $ (V2$) - (VV)' + - V) (VV)] ' dl2 + 
IF' 



+ -j^d^^}p''-'Q{±{E-V)-<!>) 



(72) 



which has been already obtained in a previous work |T2[. In practical calculations 
it is desirable to lower the number of derivatives acting on the step function. This 
can be achieved by using the identity 



V, [/(p) e(±(F -V)- $)] = 

= 2 [$ (V.$) + {E-V) (V.V)] 9^,2 [f{p) Q{±{E -V)- $)] 



(73) 



together with integration by parts yielding 
T{D) 1 



KiE) 



2 {A7ih^)D/^T{^ + 1) 



y rf^xp^e(±(F-i^) -$) 
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nHe) 
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1 



$2 



+ ^ (VV) 9^2 + -[(V.<l>) (V.V^<l>) - (v.y) (V.VV)]a^. - 



(74) 



Again, the final result for N{E) is explicitly convergent if the derivatives with respect 
to the scalar field are evaluated after space integration, as can be seen by introducing 
a new parameter A 



d^'xglx, $] 9^. [/[<l>2] e(±(^ -V)- $)" 
= J d^x g[x, $] /[<l>2 + A2] e(±(^ -V)- VWT^] 



A=0 



(75) 



It is interesting to notice that the non-relativistic limit of the former expressions 
yields that one given in section (|2.3| ). The only difference stems from the spin 
interaction with the gauge field. For instance for D = 3 the correct non-relativistic 
limit of the Dirac equation is not the Schrodinger equation but the Pauli equation 
which contains a dipole spin interaction — a - (V A A). This has its corresponding 
effect in the magnetic polarizability terms, since the pure non-relativistic treatment 
gives 



1 



48 7r2 h"^ y/2 



m 



d X 



{E-vy/ 



and the non-relativistic limit of the Dirac equation produces 



ivr^(E) 



1 



24 7r2 



d-^x 



{E - vy/' 



(76) 



(77) 



The explicit expressions for the upper branch contribution to the various densities 
are rather lengthy and are given in tables 1, 2, and 3, where the definitions 



epiyX, /i) = /i — V{x) 
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kF{x,iJ,) = - $(x)2 

XF{x,fi) = 7^ 
kp 

have been used. A few remarks are in order. The expressions have been written 
in the simplest possible way and following the standard fashion. However, they 
appear to be explicitly divergent at the classical turning hypersurface and they would 
produce non convergent mean values. In fact, when writing those expressions, some 
purely distributive terms have been intentionally omitted. They can be recovered 
by means of the rules 



Xp 



fi-V 



k 



F 



kp''Q{fi-V-<i>) ^ -2"^^dJ [kpe{fi-V-<^)] (79) 

(notice that n is always an odd number so that the derivative is evaluated an integer 
number of times). It can be checked that then any negative power of kp produces 
finite results for the mean values of the densities although the density itself still 
remains divergent at the classical turning hypersurface. Similarly to the cumulative 
number it is convenient to replace the scalar field derivatives by introducing an 
external parameter A so that the space integration can be performed first. Obviously 
such a procedure becomes rather involved in practical numerical applications. In a 
forthcoming section we will present a trick specially suited for numerical calculations 
and valid in the case of central potentials, that avoids the use of the substitution 
rules given above, hence allowing for a literal usage of tables |l|-^. 

Finally, let us remark that the case $ = m = constant, A = 0, corresponds to 
a relativistic atomic system. For such a system the semiclassical expansion up to 
h'^ order has been performed |TB| in the spirit of density functional theory. For this 



particular case our formulas in table and ^ coincide (up to purely distributional 
terms from eq. ([75|)) with their intermediate WK results. 



3.4 Application to Relativistic Harmonic Potentials 



In previous work g |Tl| harmonic potentials in the Dirac equation have been used 
to illustrate the WK expansion up to second order in h in the total energy, and a 
comparison with the exact results as well as with the Strutinsky average |3y, ^ 



have been carried out. At this order of h many of the problems discussed above do 
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not take place, in particular the distributional character of the cumulative number. 
In this section we proceed to apply our formalism to the harmonic potentials up 
to fourth order in h and for D = 3 space dimensions along the lines of |jl2|]. More 
specifically we consider the central potentials 



V{r) = $(r) — m = - mto'^r'^ , Ai{r) = . 



^0) 



The Dirac eigenvalue problem can be solved for the bound states of positive energy. 
The energy of the n-th excited state reads 



^1) 



and has degeneracy gn = (n + l)(n + 2), hence the cumulative number of states reads 



(82) 



n=0 



Straightforward application of our formulas yields 
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^3) 



where the scaled dimensionless variable x = {fi — m)/ (2m) has been introduced. 
Note that the fourth order term is a pure distribution. Similarly the total energy 
up to fourth order reads 
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It is interesting to note that the distributive contribution in the fourth order term 
corresponds to the fourth order correction of mN{^) which is the total energy of 
A^(/i) particles and hence cancels in the expression for the binding energy. We ignore 
whether this is a general feature of the semiclassical expansion. As a side remark 
let us mention that the non-relativistic reduction of the former expressions coincides 
with the non-relativistic result of section |2.4| , corresponding to the spin degeneracy 
V = 2 and a non-relativistic potential m, provided the rest mass contributions 

for the energy and the chemical potential are subtracted. 

The semiclassical expansion for this particular example can be checked directly 



as suggested in P6| , |20|| ) by means of the asymptotic Euler-Maclaurin summation 



formula 31 



n=0 



(2n)! 

(85) 

where Bm stand for the Bernoulli numbers. 

To have a quantitative idea of how large are the corrections we take m = 939 
MeV and fix the oscillator constant through the relation huj = 41y4~3 MeV. We 
also consider two nucleon species so that A = N + Z and N = Z. Our results 
are presented in table ^. We consider two versions of the semiclassical expansion 
and compare them to the exact result and the Strutinsky average. In any case we 
consider the total number of particles to be A = N{fi). In the first version of the 
semiclassical expansion, denoted as non-perturbative (NP), we proceed as follows. 
For a given order in h, say /i^^ we fix numerically the chemical potential /i directly 
from the normalization equation A = A^o(/^) + - • ■ + N2k{fJ') and recompute the energy 
up to that order as E{fi) = Eo{fi) + . . . + _E2fc(yu). This is done independently for each 
order. In the second method, denoted as perturbative (P), we proceed along the lines 
described in section (2.5) [26, 25]. We fix the zeroth order condition A = A^o(/^o) 



and compute the corrections in h to it by expanding // = yUo + + • • •• This second 
procedure has the practical advantage of making the contributions of different orders 
to be additive. We also construct the Strutinsky averaged (SA) energy since it is 
commonly believed that the semiclassical expansion should converge to it, and not to 
the exact energy since shell effects have been smoothed out. To compute this average 
we have chosen a Gaussian weight function. As can be seen from the table the 
semiclassical expansion seems indeed to converge to the Strutinsky average in this 
particular example at least up to order WK^^. We also observe a clear trend between 
the non-perturbative and the perturbative determination of the chemical potential, 
namely there appears not to be a substantial improvement in the numerical value 
of the energy. However, one has to say that the perturbative method avoids big 
cancellations as it turns out to be the case in the non-perturbative procedure, and 
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hence it is better suited for numerical estimates. Incidentally, let us mention that 
our non-perturbative calculation up to WK/i^ coincides numerically with previous 
estimates |Tl| (we only find a small deviation from our Strutinsky averag^ to theirs 
in the case A = 16). Finally, let us point out that the fourth order contribution is 
much smaller (~ 0.5 MeV) than the difference between the semi classical expansion 
and the exact result (~ 10 — 20 MeV). This is however not the case for the second 
order in h which turns out to be even four times larger than the shell effects. 

3.5 Semiclassical Dirac Sea Corrections 

As already mentioned, the Dirac sea contribution to the energy, defined as the 
sum of all Dirac eigenvalues below zero energy, is an ultraviolet divergent quantity 
which cannot be made finite by simple subtraction of the vacuum energy. However, 
the finite contributions of the sea energy acquire a precise meaning within a given 
renormalization scheme in a renormalizable quantum field theory, the infinite parts 
corresponding to absorb the infinities in the parameters of the Lagrangian. This is 
the case for instance in the a-uj model. Nevertheless, it has been shown that for 
relativistic Lagrangians of fermions interacting through Yukawa like couplings, the 
renormalization program cannot be carried out consistently within the approxima- 
tion scheme developed so far. Indeed, a short distance vacuum instability of Landau 
type takes place at the zero boson loop one fermion loop level ||3^, so that 
tachyonic poles appear. It is not clear whether this is a true feature of this kind of 
Lagrangians or simply an artifact of the approximate method used to derive it. This 
seems to suggest that it is more sensible to stick within a valence approximation. On 
the other hand, such an approximation breaks unitarity and Lorentz covariance (see 
the discussion in and references therein). Since we are dealing with an effective 
theory it is not clear what is better, to neglect the Dirac sea or to face vacuum 
instabilities. In this section we consider for future reference the WK/i^-corrections 
to the Dirac sea densities. 

From the point of view of the semiclassical expansion, such a short distance 
instability can never be seen, due to the large distance nature of the semiclassical 
approach. In addition, the ultraviolet divergencies correspond to well defined orders 
in h, namely the orders WKhP and WKh^ present quadratic and logarithmic diver- 
gencies respectively, and higher orders turn out to be convergent. In what follows 
we will derive the WKh corrections to the sea energy from the explicit expressions 
obtained within the valence approximation. First, we note that charge conjugation 

■*We stress the fact that the number of significant digits presented in our table are compatible 
with the so called plateau condition ||20[|. 
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N_{E, V) = N+{-E, -V) (86) 



synunetry implies 



Using this equation and ( ^^ - |56| ) the relations for the baryon, scalar and total energy 
densities are 

p+(x,/i, 1/) = p'{x,-fi,-V) 

p~^{x,fi,V) = -Pe(x, -/i, -1/) (87) 
respectively. In addition, the Dirac sea baryon density is given by 

P--(X) = E E C(^)^n«(x)e(-K) (88) 
n£Q a=l 

Similar expressions hold for the scalar and total energy densities. Under the as- 
sumption that the zero energy belongs to the gap we can effectively drop the step 
function and restrict the summation to the lower branch Q_. This corresponds to 
take the limit p — oo in the lower branch densities. In this limit each Dirac sea 
density has a well defined parity under the charge conjugation operation V —V. 
The former discussion amounts to the following final equations 

p««^(a;) = -p+(a;,/i = +oo) 
pTi^) = - pt (x, fx = +oo) 

pTi^) = -pJ(a;,/i = +oo) (89) 

Direct use of tables |I], ^ and ^ yields our final result for the WK^^ contributions to 
the Dirac sea densities. The limit p — oo corresponds to fc^? ^ oo and xp ^ 1. The 
densities can be looked up at table ^ for the case in which no external gauge potential 
is present. We have checked that the Dirac sea corrections to the baryon density 
can be written as a total divergence. As a consequence the baryon number remains 
unchanged. In the particular case of constant scalar field $ = m our formulas 
correspond to the sea corrections of a relativistic atomic system already considered 
in ref. |T^. We have checked that all our results agree with theirs in this particular 



case.0 In addition, we provide the scalar density. 



4 Relativistic Nuclear Potentials 



^ Notice that our definition of pE — V p corresponds to tlie density r in that reference. 
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4.1 The a-Lu Model 



In this section we analyze the semiclassical expansion in the framework of a mean 
field approximation to the a-u relativistic Lagrangian up to fourth order in h. We 
would like to stress that this model is extremely simple as a realistic description of 
finite nuclei 0, ^ |^. However, it adjusts some global properties and reproduces the 
correct shell structure. Hence we feel it is a convenient starting point to estimate 
the accuracy of the semiclassical method. To be more specific, we consider the 
Lagrangian density to be given by (we consider h = c = 1) 

L = '^[^,{td^-g,V^')-{M-g,<f))]^iI + ^{d^<f)dy- ml + 



--F^,F^'' + -mlV,V^, (90) 



1 1 
4 ^ 2 

where 

F^, = d^V,-d,V^. (91) 

The field \1/ represents the nucleon field \I' = {^^^ , 4> stands for the scalar field 
and is the vector field, and gs are dimensionless coupling constants. This 
Lagrangian may be supplemented with counterterms since it is a renormalizable 
one. We will always work in the Hartree valence approximation (i.e. no sea ap- 
proximation), where the scalar and vector fields are considered to be classical and 
the nucleus wave function is made out of a Slater determinant of A single particle 
states. We refrain in the present work from making numerical estimates of the Dirac 
sea corrections. Actually there are no Hartree calculations for finite nuclei which 
fully include the polarization of the Dirac sea to compare with. In addition one 
has the vacuum instability problems mentioned in the previous section. The time 
independent equations defining this approximation in the case of spherical nuclei 
read 

(V2 - m2)0(r) = -gsPsir) = -g^i^^)^,, 
(V2 - ml)Vo{r) = -g,p{r) = -(7,(^70^)^,1 
[-ia ■ V + g^Vo + (3{M - ^,0)] ^„±(r) = K± V^n±(r) 

V,(r)=0, ^ = 1,2,3. (92) 

Following the standard procedure we solve this set of equations self consistently 
until convergence is achieved. 

Finally the total energy of the system which corresponds to the total nuclear 
mass can be decomposed as a sum of a fermionic and a bosonic part, as follows 

E = Ep + Eq 



26 



Ef = (2ja + 1) 



E 



B 



- / d^x 



(V0)2 + m. 



- / £x 
2 



(93) 



The semiclassical expansion will be applied to the fermionic part of the total energy. 



4.2 Fixing of the Parameters 

We adjust the parameters of the model in the Hartree valence approximation to 
nuclear matter properties following [2^, This gives 

= = 357.4, C^^=gl = 273.8, (94) 

which corresponds to a Fermi momentum kp = 1.30 fm~^ and a binding energy per 
nucleon B/A = 15.75 MeV. We also fix the nucleon and vector meson masses to 
their physical values M = 939 MeV, m„ = = 783 MeV. This leaves only one 
free parameter rris which we adjust to reproduce the charge mean squared radius of 
^"Ca (3.482 fm). In summary, the set of parameters to be considered is 

M = 939MeV, m, = 449.7MeV, g, = 9.0547, 

(95) 

= 783 MeV, = 13.7997. 

In what follows we will apply only this parameter set. We emphasize once more that 
we do not pretend to make a fully realistic description of nuclear properties, but 
rather to give an idea of how large are the WK/i^ corrections to a pure Relativistic 
Fermi Gas Model. 

To have a quantitative idea of their magnitude we might plug the self-consistent 
scalar and vector potentials into the semiclassical formulas obtained in section (§). 
As it has been already mentioned, it is necessary to compute the space integrals 
after an external parameter has been introduced. For numerical computations this 
is a rather unstable procedure, since in some cases one has to evaluate up to a 
third derivative with respect to the external parameter. In the next section we will 
propose a method to circumvent this problem, at the expense of extending the radial 
dependence of the potentials to a suitable contour in the complex plane encircling 
the turning point. Strictly speaking one might do so by solving the Hartree mean 
field equations along this path. In our view this is an extremely complicated task, 
given the fact that we only want to make a numerical estimate of the corrections. 
We propose instead an alternative procedure based in adjusting the self consistent 
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potentials to a given analytical expression which can be then trivially extended to 
the complex plane. In practice we propose Woods-Saxon potentials of the form 



U 



Uo 



to adjust the self consistent potentials Uv{r 
so we make a least square fit minimizing 



l + exp{^} 

gvV{r) and U^{r) 



a 



\ 



Ell {U{rd - U{r., 
N 



(96) 



To do 



(97) 



The result of such a fit for N = 49 equidistant points between the origin and R = 12 
fm is given in table To have a more quantitative idea one can look at some nuclear 
properties such as binding energies^ and radii as presented in table |^ |^ |3^, ^ 
As we see the fit is satisfactory for the computed properties. 



4.3 Extension of the Potentials to the Complex Plane 

For the calculation of the total energy we always use the perturbative scheme, namely 
we adjust the chemical potential at zeroth order and expand around this point as 
explained in section p.5| . In general since we deal with two nucleon species we 
consider a proton and a neutron chemical potential. Furthermore, we also consider 
the radial dependence of the potentials. The basic integral appearing in the resulting 
expressions reads 

/(A) = J dr ^($, V, r) ^{fi - y (r))2 - $(r)2 - A^ 6 (^/i - V{r) - ^^rf + X^^ . 

(98) 

Other integrals can be obtained from it through successive derivatives with respect 
to the external parameter A^ at the particular value A = 0. The potentials we are 
dealing with, present a turning point at a finite value of the radial coordinate t'c(A), 
given by the equation 

/i - V(r,(A)) = v'$(re(A))2 + A2. (99) 

We assume, as it is often the case, that the function g is analytical in a region of 
the complex plane which includes the integration interval [0,rc(A)]. We also choose 
the branch cut of the square root function along the real positive axis 

1^1^ e<5-s-+-), argze[0,27r[. (100) 
^In no case have we extracted the center of mass motion. 
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With this choice the integral (0) can be transformed into an integral along the path 
defined by the lines C+ U C_ depicted in figure |I|. Under the assumption that no 
further singularities occur within the closed contour C+ U C2 U C_ U Ci it is clear 
that the original integral can be evaluated along the external line Ci of the contour, 
i.e. 

J(A) = \f dzg{<!>,V,z) ^{f,-V{z)y-^zy-X^ 

= -I dz V, z) ^{fi - Viz)y - ^zy - A2 . (101) 

In this representation, the turning points never lie on the integration path, and 
hence the integration and the derivatives with respect to the external parameter A 
at A = can be interchanged. In this case no explicit distributive contribution ap- 
pears, provided we keep away from the turning points. It should be mentioned that 
for quantities which involve logarithms there might appear in principle a different 
analytical structure. A detailed analysis reveals that also in this case formula (|101|) 
remains valid. The discussion above, allows us to readily use the results of tables 0- 
^ in the radial case literally with the additional modification of removing the step 
function and understanding the radial coordinate to lie along path Ci in the complex 
plane. Finally, it is important to check that no additional singularities appear when 
the path Ci is deformed. We do this in practice by comparing our prescription (|1U1D 



with the integral calculated along the real axis, when no derivative with respect to 
the external parameter appears. From the point of view of the complex plane, this 
condition is sufficient since derivatives in the external parameter do not introduce 
new singular points. 



4.4 Numerical Results 

We consider the Woods-Saxon potentials described above and apply the contour 
prescription explained in the previous section. The results are presented in table ^. 
There we show the numerical values for the total fermionic energy at any given order, 
together with the exact results (WS) for those potentials as well as the Strutinsky 
averaged energy (SA). The latter is computed in the cases where it can be easily 
extracted, i.e. when there are enough bound states over the chemical potential. In 
no case have we considered center of mass corrections. In general we observe sim- 
ilar trends as in the harmonic potentials case, with the exception that for lighter 
nuclei the convergence of the semiclassical series seems to be slower. It is not clear 
whether this behaviour corresponds to that of an asymptotic series. Finally, it is 
worth mentioning that in the cases where the Strutinsky average is well defined. 
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the semiclassical expansion seems to converge to a value compatible with the uncer- 
tainty inherent in the Strutinsky average. The conclusions drawn for the harmonic 
potentials seem to hold also in the case of Woods-Saxon potentials; whereas the 
WK^^ corrections to the Fermi Gas Model provide a substantial improvement to 
the total energy, this is no the case for the WK/i^ corrections, which turn out to be 
much smaller than pure shell effects. 

5 Conclusions and Perspectives 

In the present paper we have studied the role of h'^ corrections to the extended WK 
approximation within the relativistic approach to nuclear physics. For this purpose, 
we have considered the cumulative number and level density of the Dirac operator 
with scalar and vector potentials as the basic objects to which the semiclassical ap- 
proximation can be applied. This allows to compute any one-body operator within 
a semiclassical treatment and more important, reduces the calculational effort con- 
siderably. We have found extremely convenient to relate the cumulative number 
of states to the phase of the determinant of the Dirac operator. Furthermore, we 
transform the whole Dirac problem into a Schrodinger type one with some additional 
spinor structure, which makes the even nature of the h expansion more evident. In 
fact, our methodology reproduces previous calculations || rather quickly. We 
have then obtained explicit expressions for the total energy, scalar and particle den- 
sities up to ^^-order. We have found that, for central potentials, a precise meaning 
to the densities can be given by considering them along a path in the complex ra- 
dial coordinate which starts and finishes at r = and encloses the classical turning 
points. The mean values of such densities produce convergent results, at the expense 
of extending the scalar and vector potentials to the complex r-plane. We have also 
shown how the WK^^ corrections for the Dirac sea can be deduced from the valence 
densities and their explicit form has been given. 

To have a quantitative estimate of the WKh corrections, we have considered the 
relativistic a-u Lagrangian in the Hartree valence approximation. We have applied 
a semiclassical expansion to the fermionic part of the total nuclear mass and de- 
termined the contribution of the different orders in h. As external profiles we have 
considered a Woods-Saxon fit to the exact self-consistent Hartree solution, allowing 
the aforementioned extension to the complex plane. We have compared the semi- 
classical expansion with the exact result and also with the corresponding Strutinsky 
average in the cases where the latter can be easily extracted. We have found that in 
those cases the Strutinsky average produces, within its inherent uncertainty, quite 
accurately the semiclassical WK results up to fourth order. In general we have 
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found that pure shell effects are far more important than the /^^-corrections, but 
of comparable size as the contributions. Nevertheless, one should say that our 
statement is strictly true for the total energy which makes out the nuclear mass, 
but does not necessarily imply that higher h corrections are completely irrelevant 
for other observables, like e.g. form factors at not too low momentum transfers. 
It should be also mentioned that our results do not imply either that the h"^ order 
correction to the extended TF approximation has to be small. In fact, as found in 
ref. in a non-relativistic framework, the WK expansion seems to converge faster 
than a density functional approach. 

Although our method has been applied to the a-u Lagrangian, it can be directly 
generalized to more realistic and sophisticated relativistic Lagrangians including 
Coulomb interaction, p exchange, etc. which correctly reproduce a wealth of nu- 
clear data. We do not expect our results to be substantially modified by such 
generalization. 

Our numerical calculation requires further approximations besides a pure WK 
expansion. The need to produce convergent results requires an analytical continu- 
ation of the potentials to the complex radial coordinate plane. Given the fact that 
the self-consistent potentials are available on a numerical grid, we have chosen to 
fit analytical Woods-Saxon expressions to them, allowing a direct extension to the 
complex plane. One could perform better, by solving the Hartree equations on a 
given path of the complex r-plane. 

In no case have we considered a direct minimization of the semiclassical energy 
as a functional of the bosonic fields. Besides our wish to keep the problem man- 
ageable, we do not think that this is the correct way of proceeding, specially when 
the h'^ corrections are included. These corrections include fourth order derivatives 
which very strongly influence the behaviour of the system at short distances, require 
additional boundary conditions not present in the original formulation of the prob- 
lem and might even spoil the mere existence of solutions. At this point there is a 
dramatic difference with the corrections. In this case the number and form of the 
derivatives coincide with those of the pure classical bosonic energy. A systematic 
treatment of this question for higher order corrections might be achieved follow- 
ing the suggestion of ref. [25, 26 1 where the solutions are expanded in powers of h 
around the Thomas-Fermi solution, combined with our complex plane prescription. 
A detailed study of this problem is left for future research. 

In summary, we might conclude by saying that our results suggest a clear and 
compelling calculational scheme in the semiclassical approach to relativistic nuclear 
physics. The WK/i" approximation already explains many of the gross features and 
substantial improvement can be achieved by computing the next to leading order 
corrections, i.e. in the WK expansion. From there on, the semiclassical expansion 
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is not only difficult to handle, but also accounts for a tiny fraction of the h plus 
the h"^ contributions. The difference between the semiclassical expansion and the 
exact result is due to pure shell effects which are comparable in magnitude to the 
contributions and cannot be incorporated within a semiclassical treatment. 
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Appendix. Semiclassical Expansion of the Func- 
tional Determinant 



In this appendix we briefly review a very powerful method to compute a deriva- 
tive expansion of the functional determinant, which has been first suggested within 



Quantum Field Theory to compute effective actions in the low energy limit. It 
only requires the operator inside the determinant to be a differential operator of 
elliptic type, i.e. of the form 

+ U{X) (102) 

where = iTtd^ + A^{X) is the generalized momentum operator containing in the 
most general case a non-abelian gauge field An. For completeness we just quote 



the main steps of the method presented in ||2^, U{X) is an operator depending 
on the position operator and can have any internal symmetry structure. The 
following general relations are satisfied 

[X^,X,] = 
[P^,,X^] = ih6f,u 

[P^,U] = iV^U (103) 

The Euclidean indices and v run from 1 to D (the Euclidean space dimension). 
The basic object to be expanded is the determinant of the elliptic operator 

W = logDet(p2 + f/(x)) 

= Tr log(p2 + f/(X)). (104) 

In general this is an ultraviolet divergent quantity which in principle requires reg- 
ularization. However, in the applications we will be considering only the phase of 
the determinant which turns out to be convergent, and hence no explicit ultraviolet 
regularization will be introduced. In addition, the determinant is invariant under 
similarity transformations. A particularly interesting one is given by the translation 
operator in momentum space, which generates the identity 

W = Ti (es'^^ log(p2 + [/(X))e-^'-^) = Tr log ((P^ + k^f + U{X)) , (105) 

where k is an arbitrary c-numerical vector. Since this is valid for any momentum k, 
one can average over all possible values, as follows 

d^k 
d^k 



W = n-'Jj^Trlog{A-' + 2kP + P^ 

^ =lj0y^^ A(X) = (A:^ + f/(X))-\ (106) 
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It should be mentioned that the averaging procedure corresponds to project onto the 
gauge invariant subspace, and hence imphcitly assumes a gauge invariant regulator 
for W. This expression suggests an expansion in the generalized momentum operator 
by use of the formula 

W = Q-'J j^Tr log(A-i)-5:i-^(A(2A:P + p2)r . (107) 



n=l 



Although the starting average defines a gauge invariant object, the former ex- 
pansion does not preserve gauge invariance explicitly, i.e. term by term. The main 
goal is to restore this invariance in an explicit manner. We will not dwell upon the 



details which have been described at length elsewhere p8[. The final result can be 
cast in the form 



W = J -J^^ E W^2n(x, k) . (108) 



n=0 

For later reference we reproduce in table ^ the general results up to sixth order. 
Note that the infinite volume factor Q cancels in the final result, as it should be. 



34 



References 

[1] J.D. Walecka, Ann. Phys. (N.Y.) 83 (1974) 491. 

[2] B.D. Scrot and J.D. Walecka, "The Relativistic Nuclear Many-Body Problem", 
Plenum Press, New York, Advances in Nuclear Physics V.16, 1986. 

[3] P.G. Reinhard, Rep. Prog. Phys. 52 (1989) 439. 

[4] B.D. Serot, Rep. Prog. Phys. 55 (1992) 1855. 

[5] M. Centelles, X. Vines and M. Barranco, Nucl. Phys. A519 (1990) 73c. 

[6] M. K. Weigel, S. Haddad and F. Weber, J. Phys. G17 (1991) 619. 

[7] C. Speicher, R.M. Dreizler and E. Engel, Ann. of Phys. 213 (1992) 312. 

[8] D. Von-EifF, S. Haddad and M. K. Weigel, Phys. Rev. C46 (1992) 230. 

[9] M. Centelles, X. Vinas, M. Barranco, S. Marcos and R. J. Lombard, Nucl. Phys. 
A537 (1992) 486. 

[10] M. Centelles and X. Vihas, Nucl. Phys. A563 (1993) 173. 

[11] M. Centelles, X. Vinas and M. Barranco, Ann. Phys. (N. Y.) 221 (1993) 165. 

[12] E. Ruiz Arriola and L.L. Salcedo, Mod. Phys. Lett. A8 (1993) 2061. 

[13] M. Centelles in Density Functional Theory, eds. E. K. U. Gross and R. M. 
Dreizler, Plenum Press (1994). 

[14] P. Ring and P. Schuck, "The Nuclear Many-Body Problem", Springer- Verlag 
New York Inc., New York, Texts and Monographs in Physics, 1980. 

[15] M. Brack, C. Guet and H.B. Hakansson, Phys. Rep. 123 (1985) 275. 

[16] W.F. Pohner and M. Dreizler, Phys. Rev. A44 (1991) 7165. 

[17] C.H. Hodges, Can. J. Phys. 51 (1973) 1428. 

[18] D.R. Murphy, Phys. Rev. A24 (1981) 1682. 

[19] P. Sharan, Comp. Phys. Comm. 69 (1992) 235. 

[20] M. Brack and H.C. Pauh, Nucl. Phys. A207 (1973) 401. 



35 



[21] E. P. Wigner, Phys. Rev. 40 (1932) 749; for a recent review see e.g. M. Hillery, 
R. F. O'Connel, M. O. Scully and E. P. Wigner, Phys. Reports 106 (1984) 121. 

[22] E.P. Wigner, Phys. Rev. 46 (1934) 1002. 

[23] J.G. Kirkwood, Phys. Rev. 44 (1933) 31. 

[24] C.J. Horowitz and B.D. Serot, Nucl. Phys. A368 (1981) 503. 

[25] P. Schuck and X. Viiias, Phys. Lett. B302 (1993) 1. 

[26] H. Krivine, M. Casas and J. Martorell, Anna. Phys. 200 (1990) 304. 

[27] L.H. Chan, Phys. Rev. Lett. 57 (1986) 1199. 

[28] J. Caro and L.L. Salcedo, Phys. Lett. B309 (1993) 359. 

[29] R.G. Parr and W. Yang, "Density-Functional Theory of Atoms and Molecules" , 
Oxford University Press, New York, 1989. 

[30] V.M. Strutinsky, Nucl. Phys. A95 (1967) 420; A122 (1968) 1. 

[31] M. Abramowitz and LA. Stegun, "Handbook of Mathematical Functions", 
Dover Pubhcations, Inc., New York, 1970. 

[32] R.J. Perry Phys. Lett. B182 (1986) 269. 

[33] T.D. Cohen, M.K. Banerjee and C.Y. Ren, Phys. Rev. C36 (1987) 1653. 

[34] C.W. Jager, H. de Vries and C. de Vries, Atomic Data and Nuclear Data Tables 
14 (1974) 479. 

[35] D. Vautherin and D.M. Brink, Phys. Rev. C5 (1972) 626. 

[36] A.H. Wapstra, G. Audi and R. Hoekstra, Atomic Data and Nucl. Data Tables 
39 (1988) 281. 

[37] M. Centelles, M. Pi, X. Vinas, F. Garcias and M. Barranco, Nucl. Phys. A510 
(1990) 397. 



36 



TABLE CAPTIONS 



1. Wigner-Kirkwood expansion of the valence contribution to the baryonic den- 
sity, (?/'''" (a;)^/'(x))vai, in the relativistic case up to fourth order in h ior D = 3 
dimensions. The definitions ep = — V, kp = {ep — <l>^)^/^, xp = ep/kp and 
Ip = log((eir + kp)/^) have been used. For more details sea main text and 
eq. m- 

2. The same as table |l| but for the scalar density, {'ip{x)ip{x))^s,\. 

3. The same as table |l| but for the valence contribution to the total fermionic 
energy density, {'ilj\x){—ia ■ V + /5$(a;) + ^(a;))^/'(x))vai- 

4. Wigner-Kirkwood expansion for the total energy (in MeV) of a system of 
A nucleons in relativistic harmonic potentials V = ^ — m = ^mu'^r'^ where 
hu = 4:1A~^^^. WKn includes the semiclassical corrections up to /l". NP repre- 
sents the non-perturbative determination of the chemical potential, whereas P 
stands for the perturbative determination (see main text). SA is the Strutinsky 
averaged total energy with Gaussian weight factors. Eex is the exact quantum- 
mechanical total energy. In all cases the energy mA has been subtracted. The 
numbers for A, {A = Z + N and Z = N) represent the corresponding magic 
numbers for these potentials. 

5. WK^'^-corrections to the Dirac sea baryonic, scalar and total energy densities 
respectively. 

6. Numerical fit of the parameters of the Woods-Saxon Potentials (see eq. (P^)) 
to the self-consistent mean field 0, Vq, solutions in the Hartree valence ap- 
proximation with 49 points, a (in MeV) is the standard deviation between the 
Hartree solutions and the Woods-Saxon potentials within the physically rele- 
vant region < r < 12 fm. Uo is given in MeV, and Ro and Oo are presented 
in fm. 

7. Comparison of the binding energy per nucleon, B/A (MeV), and mean squared 
charge radius, m.s.c.r. (fm), for several closed shell nuclei, a-u means the a-LO 
self consistent Hartree solutions with the parameter set specified in the text. 
WS stands for the best fit of Woods-Saxon potentials to the corresponding self 
consistent scalar and vector fields. For general information we also quote the 
experimental values whenever they are accessible. 

8. Semiclassical expansion for the total fermionic energy in MeV for several closed 
shell nuclei. The scalar an vector potentials are of Woods-Saxon type fitted 
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to the self consistent Hartree valence potentials (see table |]). The chemi- 
cal potential has been always adjusted within the perturbative method (see 
section p.5|) . WK„ contains the corrections up to order h'^. SA stands for 
the Strutinsky averaged energy and WS represents the quantum mechanical 
energy. In all cases the rest mass of the constituent nucleons M has been 
subtracted. 



9. Values of the coefficients W2n{x,k) as defined in eq. ( |108|) . The following 
convention has been used: 

(a) Given an object Yj with / = (ai . . . a„) a set of Euclidean indices we 
define Yj^ to mean —i[P^,Yj]; for instance A^,^ = {—iY[P^,[P^, A]] ; 

(b) Each term in the table means the half-sum of the term itself and its 
mirror symmetric. The mirror symmetry is defined as —>■ P^ = P^; 
A ^ = A and ^ (YZf = ZY. For instance A J = -A^; 

(c) Notice the cyclic permutation symmetry. 



FIGURE CAPTION 

1. Contour in the radial complex plane used to evaluate radial integrals in the 
classically allowed region (defined by the path C+), namely J(j^ = ~| /cr The 
branch cut starts at the radial turning point rc(A) defined in the text. 
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+ ^ (3 - 


^ (2 - xf') {y^f + Ixp (3 - X, 
')(yvf - ^<f(V2$) - 2{xF + 2lF)iy^V) 


.2)(V.$)(V.r) 

>e(/i-F-$) 



P4=2^ T7ii^'^'(V''<f) + ^xH-5+x^2) (W) 



+^x^(15-10xi.2 + 3xi.4) (V.V2V)(V.$) 
i-^<l> (-14 + 15x^2) (v,$)(V,V,$)(V,<l') 



30 

+^$=^(V,V^)(V,V,y)(V,-$) + ^ (96 - 280x^2^ 175 x/) (V-l-)^ 
U (900 - 2455 xf^ + 3351 xp^ - 2145 + 525 x/) ( Vi$) ( Vil/) ( V$)^ 



24 TT^ j lOfcj?'^' V ^ ' 10 # 

^(-5 + 4x^2) (v,V,<l>)' + -1-^ (-19 + 12x^2) (V.V.F) 

■^xf (45 - 35x^2 ^ ^2x/) (ViV,$)(ViV,r) 

3 ^ /'c c ^ 2 



10 



XF (5 - 5 xf^ + 2 x/) (Vil^)(ViV2$) 



L (-7 + 6x^2) (v,$)(V.V2$) + 1 (-5 + 3x^2) (Vil^)(V.W) 



10 fcF'' 

^ xf (-75 + 125xi.2 _ 102 x/ + 30 x/) (Viy)(ViVj<l>)(Vj$) 



^Xir(-135 + 205xF^ - 159xf^ + 45xf^) (Vi$)(ViVjV)(Vj$) 



(18 - 55x^2 ^ 35 ^^4) (y^)2^^^^2 
Xf (-54 + 73 Xi.2 - 54 x^^ + 15 xf^) ( W) ( V$)^ 



12 



1-^$ (-13 + 15x^2) (v,\/)(V,V,$)(V,l^) 



lOA.'F'' 

^ xj. (-90 + 185 xf^ - 156 xf^ + 45 x/) ( V^l/) ( V^V ) ( V,l^) 
^ (89 - 275xi.2 + 175x/) {V ,V){V ,V) 
-^xf (675 - 2270 xf^ + 3264 x^^ - 2130 x/ + 525 x/) (VV)^(Vi$)(Vi1/) 

I ^ / o« , OK ™ 2N /v72i 



^<l>(-26 + 25xp^) (W)(Vi$)(Viy) 



10 A:f 

^ (63 - 270 xp2 + 3-/) (vy)^ 
^$(-18 + 25x^2) (V$)2(V2$) 

^xj.(-135 + 265x^2 _ 237a;/ + 75a;/) (Vi$)(ViF)(V2$) 

1 rf, /' Q I r:^_2\ /'V7T/^2/V72/IS^ i 1 /' Q i c;^_2 



^,^-,$(-3 + 5x^2) iyVYiy'^) + ^(-3 + 5x^2) (V2<l>)2 

+ ^a;i.(9-8xf2^3x/) (V2r)(V2$) 
^ Xi. (-135 + 305 x^2 _ 261 x/ + 75 x^^) (VV)^(V2V) 

^(-26 + 15x^^)(VV)^ + ^4 |0(/.-r-$) 
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Ps 2 24 7r2a 


(2 + (V<l>)^ - ^ <I>(V2\/) - ^ (2 + X, 
=f') (Vy)' + 2(-a;^ + 3/^) (V^^) 


^2) (V.<l>)(V,r) 

.e(/i-y-$) 



p^4 = 2^|toW^^(2+v)(v^<^') + Td^^'^i'^'y) 

1 n I /I ^_2\ /V7 V7 ^^^V7 V7 T/^ I 1 ^ I OK ^_2\ /V7T/'\2 



+ 5|3Xf (-2 - Xi.^ + (v,$)(V.V2<|.) 

1 I Q ™ 2\ ^V7 TA\^V7 V72/R\ i 1 



(l + AxF^) (V.V,$)(V.V,V) + ^$(3 + 25x^2) (v\/)2(v2v) 



1 ^ / n 1 2 I 1 o ^ 4\ /V7 V7 T^^2 , 1 



xf {-9- 17 xf^ + 12 xf^) {ViVjVy + ^ x^.^ (VV)(V2$) 



^(2 + 3xi.2) (V.V^)(V.V2$) + (-1 + 6x^2) (v,VV)(V,<l') 



^x^ (-1-3x^2^2x^.4) (V.V)(V.W) 

XiT (60 + 35 XF^ - 96 x/ + 45 x/) ( V,$) ( Vi V, $) ( $) 
^$(2 + 5x^2) (V,V^)(V,V,$) (V,*) 

^$(4+15X^2) ^^^^^^^^^^y^^^^^^ 



^Xi.(15 + 35xir2 -66xf4 + 30x/) (Viy)(ViV,l^)(Vj$) 
■gQ^ xf (-160 - 80 x^^ + 507xir4 - 530 x^^ + 175 x/) (V$)^ 
(_8 _ 18x^2 ^ 35 ^^4) (v,$)(ViV^)(V$)' 

^ Xi. (-30 - 85 xf^ + 309 x^^ - 315 xf^ + 105 x^.^) {VVf{V<^f 



24*5 

XF (30 + 40 XF^ - 93 XF^ + 45 x/) (ViV) ( Vi V, $) ( V, 
^$(1 + 3x^2) ^^^^^^^^^^.^^(^^.^^ 



^ Xf (-165 - 415 xf^ + 1542 x^^ - 1575 x/ + 525 xf®) ((Vi$) (Vi^))' 
(-6 - 17x^2 + 35x^4) {VVf{Vj<l>){VjV) 



■^xf^ (130 - 183 Xf' + 75 x/) {V^V)iVj^){VjV) 

Xf (-135 - 400 Xf' + 1506 xf^ - 1560 xf*' + 525 xf^) (Vl^)^ 
Xf (120 + 50 Xf' - 159 xf^ + 75 xf^) (V$)' (V'<I>) 
^ $ (11 + 25 xf') (V,*) (V.\/) (V^) 

Y2^xf (6 + 13xf'-30xf4 + 15xf^) (Vy)'(V'$) 
^XF (-33-14xf' + 15xf') (V'<I')(V'<I>) 

gO^XF(18-31XF' + 15XF^) (V'V)' + ^XFF^\e{i2-V-^) 
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pl, = ^{elkF + eFkl-^Hp}Q{^i-V-^)+ptV 

= 24^|2 $ {-xp + h) (V2$) - 2 ^ (V.$)(V.r) - 2 A;^ (1 + x^^) {V'V) 
+ {xf-xf^ -lF){y^f + {2xf-xf^ -2lF){yvf^Q{p-V + ptV 

Pe. = ^1 Th^FHV^<^) + j^{-3 + xf'){W) 



10 kp 

'+12 

1 i'_o _L /I ^„2WY7.Y7 .T/\2 ' 1 ^_3 



(_23 + 12 (V, V,*)' + $ (-7 + 12 (y, v,<l>) (V. V,\/) 



20*2^^ (-9 + 4 V) (V.V.y)^ + j^xfH-II + Qxf') (V,<l>)(ViV2$) 



10 feF"* ^ ^ ^ * ^ ^ * ^ 5 /CF 

-^^Lxi.^ (100 - 123x^2 ^ 45 ^^4) ^^.^^^^.^^.^^^^.^^ 
(1 - 36x^2 + 30 3,^4) (v,y)(V,V,<l')(V,$) 
(1 - 12x^2 + 93;/) (v,$)(V,V,y) (V,$) 



ekp'^ 

x^^ (_208 + 444 Xi.2 _ 363 x/ + 105 x/) (V$)^(V$)^ 



+^xf' (35 - 42x^2 ^ 15 ^^4) (v,\/)(V,V,\/)(V,<l') 



$ (-6 - 165 x^2 _^ 175 XF^) ( V.<f ) ( V.F) ( V<l>)^ 
■^XF^ (-40 + 87xi.2-72xf4 + 21x/) (V\/)'(V<I>)' 
^ (4 - 21 x^2 + 15 a;^4) (vV)(V<l>)^ 



12A;f 

+^ x^3 (19 _ 24x^2 + 9 3-^4) (v,l^)(V,V,$)(V,F) 



(_8 - 57X^2 ^ 45 ^^4) ^^^y^^^^^^y^^^^y^ 

xf^ (-995 + 2172 xf^ - 1800 xf^ + 525 x/) (( V,$) (Vil^))^ 
$(_19_ 160x^2^ 175 x/) (V\/)'(V,$)(V,V) 
x^^ - 71 xf' + 25 x/) (W)(V,$)(V,\/) 
Xf 3 (-285 + 689 xf^ - 595 x/ + 175 x/) ( VV")^ 
xf^ (140 - 189 xf^ + 75 x/) (V$)^(V2$) 
K (-1 - 24xi.2 + 25xi.4) (V.$)(Vi^)(V2<l') 



10 A;f 

^ x^3 (25 _ 36 XF^ + 15 x/) (V\/)'(V2$) 

1 ^_3 /'_10 _L 1 ^_2WV72dS^2 , 1 



^Xf3(_19 + 15 3,^2) (y2$)2 ^ 1 $ (_1 + 3 a;^2) (y2y) (y2$) 



60$^ V y V " -7 ' 6fcF'=' 

(-19 - 96x^2 ^ 75 ^^4) (V\/)2(V2\/) 

12 + 5xi.2) (V2V)2 + ^(a;^-/^)Fi/ |e(;U-V-<l>) + p|l^ 



20* 



1 s 
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A 


WKo WK2 WK4 


SA 


Ecx 


4 


NP 
P 


139.428 162.018 161.104 
139.428 160.405 161.066 


161.1 


153.412 


16 


NP 
P 


557.711 592.013 591.437 
557.711 591.011 591.427 


591.44 


579.979 


40 


NP 
P 


1394.28 1440.21 1439.78 
1394.28 1439.47 1439.78 


1439.78 


1424.52 


80 


NP 
P 


2788.56 2846.08 2845.74 
2788.56 2845.50 2845.74 


2845.74 


2826.66 


140 


NP 
P 


4879.97 4949.07 4948.79 
4879.97 4948.59 4848.79 


4948.80 


4925.89 


224 


NP 
P 


7807.96 7888.63 7888.39 
7807.96 7888.21 7888.39 


7888.39 


7861.67 



Table | 
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-f|$-^(v.$)(v,i^)(v$)^ + |$-^(W)(v$)' + ^$-4(v,y)(v.v,i^)(v,\/) 



+|<I>-3(V,<I>)(V,V2$) + |$-3(V,y)(V,W) - f|$-4(V,<l>)(V,V,<f)(V,<l>) 
-|$-4(ViV)(V,Vjl^)(Vj$) + ii$-5(v<l>)^ + |$-5(VV)^(V$)^ 
-ii$-^(V.r)(V.V,<l>)(V,\/) + f|$-M(V.$)(V.\/))^ - il<|.-^(V2V)(V,$)(V,V) 
+^$-5(vi/)4 - f $-4(V$)'(V2<l>) - i$-4(VV)'(V2$) 



+ i<l>-2(V.$)(ViV2$) + f $-2(y.y)(y.y2y) _ li $-3 ( V,<|.) (V, V, $) ( $) 

-|<l>-3(Viy)(V,Vjy)(Vj4>) + ii<l>-4(V$)^ + |$-4(V\/)^(V$)^ 
-I <l>-3(V,F)(V.V,<l>)(V,y) + f ((V.$)(V.\/))' 
_Z$-3(vV)(Vi$)(Vil^) + i<l>-4(Vy)' - |§$-3(V<I>)'(V2<I>) 

-|<I>-3(VV)2(V2$) + ^$-2(V2$)2 + X$-2(V2V)2| + psea^ 



Table § 
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Ro 






a 




V 


343.8 


2.630 





5602 


2.33 






415.4 


2.624 





6156 


2.16 


20 Ca 


V 


401.5 


3.518 





6831 


5.76 







478.8 


3.550 





7175 


5.19 


20 Ca 


V 


386.7 


3.880 





6074 


5.35 







464.5 


3.902 





6458 


4.52 


iNi 


V 


379.2 


4.175 





5265 


4.97 







457.9 


4.190 





5698 


4.57 




V 


374.0 


4.954 





5752 


5.28 







452.2 


4.973 





6193 


5.21 




\/ 


371.7 


5.685 





5745 


1.86 







449.6 


5.712 





6151 


2.88 




V 


377.0 


6.634 





6200 


1.99 







455.0 


6.669 





6571 


1.94 



Table ^ 
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a-u 


WS 


Exp. 




B/A 


m.s.c.r. 


B/A 


m.s.c.r. 


B/A 


m.s.c.r. 




3.33 


2.84 


3.38 


2.83 


7.98 


2.73 


40 p„ 


6.27 


3.48* 


6.22 


3.48 


8.55 


3.48 


20 


6.50 


3.47 


6.34 


3.49 


8.67 


3.47 


iNi 


7.23 


3.72 


7.00 


3.73 


8.64 






8.35 


4.22 


8.10 


4.21 


8.71 


4.27 




8.80 


4.60 


8.58 


4.63 


8.36 






9.83 


5.34 


9.59 


5.37 


7.87 


5.50 
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